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ABSTRACT 

An  existing  program,  which  determines  the  inverse  Laplace  transform  of  a 
quotient  of  two  polynomials,  provides  expanded  systems  analysis  capabiltiy.  The 
program  complements  a  Root  Locus  Program  (AMSAA  Technical  Memorandum  No. 
21)  and  a  Frequency  Response  Program  (AMSAA  Technical  Memorandum  No.  69). 
The  program  uses  a  self-contained  complex  arithmetic  routine  and  also  a 
self-adjusting  variable  scale  plotting  technique.  The  plotting  is  done  on  a  standard 
line  printer  and  gives  a  time  history  plot  of  system  response  for  a  variety  of  input 
forcing  functions. 

A  listing  of  the  FORTRAN  IV  source  deck  and  the  corresponding  flow 
chart  of  the  program  is  shown  in  the  appendixes.  Also,  several  examples  are  given 
to  introduce  the  user  to  the  operating  procedures  and  capabilities  of  the  program. 
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1.  INTRODUCTION 


One  method  of  analyzing  physical  systems  requires  the  use  of  Laplace 
transforms.  The  inverse  Laplace  transformation  is  necessary  to  analyze  a  system’s 
time  response  to  a  forcing  function. 

Time  constants,  overshoots,  resonant  frequencies,  and  steady  state  errors  are 
obtained  by  the  time  response  of  the  system.  By  observing  these  characteristics,  the 
designer  is  able  to  optimize  the  system  response  to  a  particular  performance  index. 

Inverse  Laplace  transforms  can  easily  be  found  in  a  handbook  of  Laplace 
transform  pairs  for  the  simpler  control  system  equations.  For  more  complex  systems 

it  is  a  tedious  task  to  take  the  partial  fractions  representing  the  quotient,  find  the 

corresponding  inverse  Laplace  transform  of  each  individual  term  and  sum  them. 

An  existing  FORTRAN  IV  program  (References  1,2,3)  has  been  adapted 
for  system  analysis  application  to  relieve  the  computational  burden.  The  program 
determines  the  inverse  Laplace  transform  of  a  wide  variety  of  high-order  system 

equations.  The  transient  response  is  graphically  displayed  for  quick  and  easy  analysis. 

The  main  features  of  the  modified  program  are: 

1.  FORTRAN  IV  programming  which  needs  no  machine-oriented  or  object 

language. 


2.  Self-contained  complex  arithmetic  (compatible  with  BRLESC). 

3.  Digital  plotting  with  on-line  printer. 


HlcCracken,  Daniel  D;  FORTRAN  With  Engineering  Applications;  John  Wiley 
and  Sons,  Inc.,  New  York,  1967,  pp.  196-214. 

2 

Titus,  C.  K. ;  A  General  Card  -  Program  for  the  Evaluation  of  the  Inverse 
Laplace  Transform;  J.  Assoc.,  Computing  Machinery,  2_>  No.  1,  January  1958. 

3 

T.  R.  Boshkow,;  Curve  Plotting  Routine  for  the  Inverse  Laplace 
Transform  of  Rational  Functions;  J.  Assoc.,  Computing  Machinery,  .5, 

No.  1,  January  1958. 
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4.  Maximum  order  of  numerator  and  denominator  polynomials  are  51  and 
52  respectively. 

5.  Fast  operating  time. 


Main  features  of  the  graphical  display  are: 

1.  Automatic  scaling  to  give  maximum  viewing  attenuation. 

2.  Preselected  scaling  to  give  uniformity. 

3.  Variable  self-adjusting  computing  interval  for  quick  running  time. 

4.  Manual  selection  of  up  to  15  different  real-time  intervals  for  a 
particular  function. 

5.  Automatic  cut  off  at  time  when  response  is  constrained  within  1  X 
10"^  for  five  consecutive  computation  intervals  (where  K  is  dependent  upon  the 
specific  problem  being  solved). 

Requirements  of  the  Program 

The  quantities  needed  for  problem  solution  are: 


1.  Coefficients  of  numerator  and  denominator  of  the  system  transfer 
function  in  polynomial  form^’^ 


N  (s) 

D(s) 


n  ,  n-1  n-2 

as  +  a,s  +  a0s  + 
o  1  z 


a  .  s 
n-k 


,  m+n  ,  ,  m+n-1  ,  ,  m+n-2  ,  ,  J 

b  s  +  bns  +  b0s  +  .  .  .  b  .  fls 

o  1  2  m+n-£ 


(1) 


4 

Burke,  Harold,  H.  and  Payne,  Robert  L.;  A  Linear  Closed  Loop  System 
Analysis  Procedure  Using  Line  Printer  Plots  of  Characteristic  Equation 

Root  Loci,  AMSAA  TM  No.  21,  November  1968,  Army  Materiel  System  Analysis 
Agency,  Aberdeen  Proving  Ground,  Maryland. 

^ Burke,  Harold,  H.  and  Payne,  Robert  L. ;  A  Linear  and  Nonlinear  Systems 
Analysis  Tool;  Line  Printer  Plots  of  Characteristic  Equation  Root  Loci, 

Bode  and  Popov  Plots  of  System  Transfer  Functions,  AMSAA  Technical 
Memorandum  No.  69,  Army  Materiel  Systems  Analysis  Agency,  Aberdeen 
Proving  Ground,  Maryland.  (In  publication). 
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2.  Roots  and  the  multiplicity  of  the  roots  of  the  denominator. 


4,5 


2.  THEORY 


The  basic  transfer  function  of  a  closed-loop  system  can  be  shown  to  be  of 
the  form  given  in  Figure  1.  For  multiple-loop  system,  the  G’s  and  H’s  are  readily 
expressed  as  sums  of  products  of  polynomials  which  are  identified  with  the 
individual  elements  making  up  the  complete  system. 


Kc  H |  (S) 
H2  (S) 


Figure  1.  Linear  Closed-Loop  System 


where 


R  «  System  Input 

C  =  System  Output 

F  =  System  Feedback 

E  =  System  Error 

K  =  Process  Gain 
P 

K  =  Controller  Gain 
c 


The  linear  open-loop  system  transfer  function  is 


F 

E 


G1<S>I 

KP  G2(s)J[Kc 


H^s) 

H2(s)_  • 


(2) 


4Ibid. 

5Ibid. 
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The  linear  closed-loop  system  transfer  function  is 


KpG1(s)H2(s) 


R  K  K  G  (s)H.(s)  +  G.(s)H,(s)  * 

c  p  1  1  2  2 


(3) 


Equations  (2)  and  (3)  can  be  expressed  as  ratios  of  polynomials  as  shown  in 
Equation  (1).  Rewritten  in  series  form,  Equation  (1)  becomes 


n-k 

V 

N(s)  _  jjSp 


a  s 
V 


n-y 


D(s)  m-t-n-£ 


(A) 


biS 


m+n-i 


where  the  numerator  and  denominator  polynomials  are  different  for  the  open-  and 
closed-loop  system  transfer  functions. 


The  factored  form  of  the  right-hand  side  of  Equation  (4)  is 

sk  n  (s+z  )  n  fs  +  (a  +ja)  )1  [s  +  (o  -ju  )1 
U=1  M  q=v+l  L  44JL  44J 


£ 


n  (s+p.)h  n  fs  +  (o  +jci>  )]  [s  +  (o  -ju  )] 
i=l  t-g+l 


(5) 


where  h  =  1  or  2 

n  -  k  -  v 


,  and 


m  +  n  -  £  -  q 
r  “  2 
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The  inverse  Laplace  transform,  c(t)  for  an  arbitrary  input  forcing  function  r(t)  of 
Equation  (4)  is 


c(t) 


R(s) 


k  v  y 

s  n  (s+z  )  n 

u=l  M  Q=V+1 


[s+(oq+jUq)][s+(oq-jUq)J 


n  (s+p  )h  n  f"s+(a  +ju  )"|  |"(s+o  -ju  )1 
i=l  t-g+1  L  JL  J 


(6) 


where  R(s)  is  the  Laplace  transform  of  the  input  forcing  function  r(t).  R(s)  will  be 
multiplied  by  the  right-hand  side  of  Equation  (4)  to  obtain  the  Laplace  transform 
given  in  Equation  (6).  The  partial  fraction  expansion  of  Equation  (6)  is 


c(t) 


1  v-  g 

E  VE 


K 


ih 


J-A-k 


i-1  (s+Pi> 


E 


K 


2  2 
, ,  s  +2acjs  +  u 
t=g+l  t  t  t 


,  (7) 


where  Kj,  Kj^  and  Kt  are  the  residues  at  each  denominator  (pole)  singularity  of  the 
right-hand  side  of  Equation  (6)  and  h  =  1  or  2 


The  coefficients  of  non-repeated  real  roots  are  at  origin 

Kj  ‘  Llmlt  [S  R(f}_  > 

j=l  s=o 

on  real  axis 

Kih  ‘  Llmlt  [(S+Pi)  R(f}]  • 

h-1  s=-pi 


(8) 


(9) 
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The  coefficients  of  the  partial  fraction  terms  of  Equation  (7)  related  to 


the  repeated  real  roots  are 


K  =  f_JL _ £isi\\ 

Kj  | (j-d*  daj-i  L(s  ) 


(10) 


s=o 


and 


K 


ih 


(  1  dh_1  r  h  C(s)l| 

|  (h-1) !  dgh-l  L(s+Pi)  R(s)J) 


(ID 


s=-p 


where  h  =  2. 


The  coefficients  of  the  partial  fraction  terms  of  Equation  (7)  related  to 
the  complex  roots  are 


QSsl 


\  ■ llmlt  I  tgf [s  +  (° 

s  =  “(ot+jut) 


(12) 


and 


\2  ■ llmlt  («w  [s  +  (0 

S  =  -(o  -jw  ) 


(13) 


where  Kt  in  Equation  (7)  equals  Kt^  +  Ktv 

Once  these  residues  (Kj,  and  K^)  have  been  determined,  the  time 

response  of  c(t)  can  be  determined. 


c(t) 


A-k  v  .(j-1)  R  y  (h-1) 

_  +  V  ihfc 

j=1  (j-1)!  fr{  (h-1)! 


-E  3 


e-pi 


'  +  E  \ 

t=g+l  l 


e-(a+ j«t) 


+  K  , 

Z2 


h  =  1  or  2 


(14) 
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or  the  transcendental  form  of  Equation  (14)  is 


*-k  K  .(j-l)  g  „  (h-1) 


(h-1)! 


j-1 


i=l 


iC  +  ^  2lKtk  0t  sin  («ott+^Kt)  . 
t=g+l 


h  =  1  or  2 


(15) 


3.  GRAPHICAL  METHODS 

A  plotting  routine  is  used  to  graphically  show  the  time  response  of  the 
inverse  Laplace  transform.  Equation  (15).  The  routine  varies  the  solution  interval 
in  accordance  with  the  specific  problem  being  solved.  This  adjustable  time  step  is 
computed  as  a  function  of  previous  step  size  and  the  slope  between  consecutive 
points. 


The  variable  time  step  is  initiated  by  first  finding  the  smallest  time 

increment  after  the  problem  start  time.  This  initial  step  is  a  function  of  the  largest 

real  root  (fastest  exponential  rise)  and  of  the  largest  imaginary  root  (highest 

frequency  sine  wave).  The  increment  thereafter  will  be  determined  by  the  change  of 

the  slope  between  points.  The  expression  for  the  derivative,  f(t),  is  found  by 
explicit  differentiation  of  the  function.  The  time  step  is  then  determined  by  means 
of  comparison  between  the  slope  <£j  =  tan‘*f(t)  evaluated  at  t..  and  the  slope  <{k  +  1 
evaluated  at  tj  +  Atj.  The  solution  time  is  left  to  the  discretion  of  the  programmer. 

Critical  areas  may  be  displayed  by  isolating  the  particular  region  of  interest. 

There  will  be  a  priority  cutoff  time  if  the  response  is  constant  within  a 
difference  of  1  x  10'^  for  five  consecutive  computing  intervals.  (K  is  equal  to  8  in 
the  current  program.)  The  response  will  then  be  assumed  constant  out  to  the  stop 
time  requested  in  the  input  data.  This  procedure  is  employed  to  allow  for  maximum 
sensitivity  of  the  line  printer  plotting  display. 

It  is  not  necessary  to  specify  the  upper  and  lower  limits  of  the  c(t)  axis. 

The  program  will  select  a  scale  which  permits  the  graph  to  be  displayed  with 
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maximum  sensitivity.  If  no  limits  are  specified,  the  program  will  adjust  the  c(t)  axis 
scale  dependent  upon  the  maximum  and  minimum  values  of  the  response  function 
contained  in  the  computation  interval. 

Parametric  plots  may  be  directly  compared  by  using  the  same  c(t) 
amplitude  scales.  Only  the  values  contained  in  the  region  bounded  by  the  specified 
limits  of  the  c(t)  axis  will  be  plotted. 


4.  DATA  FORMATS 

4.1  Inverse  Laplace  Transform. 


1.  General  Description. 


a.  Computes  the  inverse  Laplace  transform  parameters  of  the  first,  second 
and  third  order  poles  at  the  origin.  No  input  data  are  necessary  for  this  process 
because  it  is  evaluated  automatically. 

b.  Computes  the  inverse  Laplace  transforms  of  first-  and  second-order  real 
poles  not  at  the  origin. 

c.  Computes  the  inverse  Laplace  transforms  of  the  first-order  complex 

poles. 


2.  Input. 

Description 

Columns 

Data 

(a)  Coefficients 

1-5 

Order  of  term  of 

of  Numerator 

(Integer) 

polynomial  (up  to  51st 

Polynomial 

order) 

6-20 

Coefficients  of 

(Floating  Pt.) 

corresponding  term 

(b)  Coefficients  1-5  Order  (up  to  52nd) 

of  Denominator  (Integer) 

Polynomial 

6-20  Coefficients  of 

______________  ^Floating  Pt _  _  _corresponding_term 

16 


Description 


Column 


Data 


(c)  Roots  of  Denominator  1-15 

(Floating  Pt.) 


Real  Part  (a^  or  p^) 


and  (p±) 


Poles  of  forms 


(16-30) 

(Floating  Pt.) 


Imaginary  Part  (u^) 


(31-45) 

(Integer) 


Multiplicity  of  Root 


3.  Output. 

The  output  of  the  inverse  Laplace  transform  consists  of: 

a.  Real  part  of  root  (a^  or  pj). 

b.  Imaginary  part  of  root  (cot). 

c.  Multiplicity  of  root  (K). 

d.  Inverse  Laplace  transform  equation  of  each  individual  partial  fraction. 

4.  Restriction. 

a.  Maximum  order  of  numerator  is  51. 

b.  Maximum  order  of  denominator  is  52. 

c.  Maximum  orders  of  K,  multiplicity  of  roots  are: 

Roots  at  origin,  0  <_  K  <_  3. 

Roots  not  at  origin  on  real  axis,  0  <_  K  2. 

Roots  in  complex  plane,  0  <_  K  <_  1. 
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4.2  Plotting  Routine  for  Inverse  Laplace  Transform. 


1.  General  Description. 


a.  Computes  the  numerical  values  of  the  inverse  Laplace  transform  as  a 
function  of  time. 


b.  Produces  a  graphical  display  using  line  printer. 


2.  Input. 


Description 

Columns 

Data 

(a) 

Identification 

1-80 

(Alphanumeric) 

Identification  of 
problem 

(b) 

Running  time  of 
problem 

1-15 

(Floating  Pt.) 

Start  time 

16-30 

(Floating  Pt.) 

Stop  time 

Scaling  of 

Output  Response* 

(31-40) 

(Floating  Pt.) 

Minimum  value 

(41-50) 

(Floating  Pt.) 

Maximum 

3.  Output. 

Produces  a  time  response,  tabulated  output  and  graphical  display  using  on¬ 
line  printer  consisting  of: 

a.  Start  and  stop  time  of  computation. 

b.  Automatic  or  preselected  scaled  ordinate  axis  with  values. 


* 

The  output  response  is  automatically  scaled  if  minimum  and  maximum  values 
are  not  specified.  ,ft 


c.  Automatically  scaled  time  axis  with  values. 


d.  Response  plots  of  tabulated  output. 
4.  Restrictions. 


a.  Cannot  exceed  1 5  time  intervals  of  a  particular  problem  without  reading 
in  coefficients  again. 

b.  Independent  variable  time  must  be  positive  (trivial  requirement  of 
physical  systems). 


4.3  Placement  of  Data  Cards. 

Group  No,  Content  Description  Number  of  Cards 

1  Identification  of  Data  1 

2  Time  (start/stop  ,  Min/Max  Val)  Max  15 

3  Blank  Card 

4  Power;  Coefficients  of  Numerator  Max  51 

5  Power;  Coefficients  of  Max  52 


Power;  Coefficients  of 
Denominator 

Roots  of  Denominator, 
Multiplicity 

Blank  Card 


As  many  as  needed 


Data  placement  is  repeated  for  next  problem  making  certain  that  the  blank 
cards  are  properly  placed. 


5.  CAPABILITY  OF  THE  PROGRAM 

The  type  of  input  forcing  function  that  the  system  is  subjected  to  must 
be  included  in  the  ratio  of  polynomials  of  Equation  (4).  A  new  polynomial  ratio  is 
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formed  by  this  multiplication. 


g(s)  =  gj(s)g2(s), 

where  gj(s)  is  the  system  input  forcing  function,  and  g2(s)  is  the  system  transfer 
function. 


The  program  is  capable 

of  handling  the  following  system 

input  forcing 

functions,  gj(s). 

F(t) 

Type  of  Input 

g-j^s) 

6 

impulse 

1 

1 

step 

1/s 

t 

rate 

1/s2 

t2/2 

parabolic 

1/S3 

-at 

e 

exponential  lag 

1 

s  4*  a 

t  -t/T 

2  e 

2nd  order 

i 

TZ 

exponential  lag 

(s  +  a)2 

2  2  1/2 

a)(l+a  a)  )  sin 
\p  =  tan  ^  (aw) 


Undamped  phase  1  t  as 

shifted  sinusoid  - 

1  2 

0) 


0) 


[1  -  2a£o)  +  a^o)^ 

1  -  c2  -I 


2  2l  1/2 


„■»_  /  /l-JT 2‘ 


sin  (a)/l-£  t-hj/) 


tan 


-1  (atp/l-Q 
1  -  a£a) 


Damped  phase 

shifted 

sinusoid 


1  +  as 


1  + 


0) 
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The  inverse  Laplace  transform  digital  program  is  capable  of  solving  any 
system  transfer  function  being  acted  on  by  a  forcing  function  which  can  be  broken 
down  into  the  following  partial  fraction  forms: 

Location  of  Denominator 

_ Roots  (Poles) _  Partial  Fraction  Form 


Origin  c^  c_  c_ 

s’  2’  3 

s  s 


Real  Axis 


c  c 

s+a  *  (s+a) 


2 


Complex  Plane 


2  2 
(s  +  2au)s  +  w  ) 


The  maximum  order  of  the  numerator  of  g(s)  can  not  exceed  51  and  the 
maximum  order  of  the  denominator  of  g(s)  can  not  exceed  52. 


6.  EXAMPLE  PROBLEMS 


Three  transfer  functions,  g2(s),  are  subjected  to  various  input  forcing 
functions,  gj(s),  to  demonstrate  the  utility  of  the  program.  Tables  1,  2,  and  3 
demonstrate  the  mirror  of  the  input  for  each  system  subjected  to  one  specific  input 
forcing  function. 

®2N 

EXAMPLE  1:  g  (s)  =  -£il 

1  g2D 


where : 


g2N  =  3s  +  11 
g2D  s2  +  8s  +  15 


g^s) 


1  1  1  +  0.5s 

o  »  o  +  1  »  2 

S  X  1  +  0.3s  +  0.25s 
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EXAMPLE  2: 


where : 


g 

EXAMPLE 

where : 


*2N 


82D1  +  82D2  +  82D3  +  82D4  +  82D5  +  82D6 


®2N 

82D1 

g2D2 


1.34  x  10~3s2  +  2.067  x  10_1s  +  1.0 


6.9981  x  10~7s5 


1.15529  x  10  4sA 


g 


2D3 


’2D4 


g 


2D5 


1.68  x  10  V 


3.828  x  10  3s2 


2.593  x  10  1s 


g 


1 

3: 


2D6 

(s) 


1.0 


s’  s  +  1  *  2  ’  2 

(s+ir  1  +  0.2666s  +  O.llls^ 


,  .  g2Nl  +  g2N2  +  S2N3  +  g2N4  +  g2N5 

^2  S  g  +2  +£  +2  +  2 

°2D1  g2D2  82D3  82D4  82D5 


govn  =  2,46665  x  lO^s^ 
2N1 

g2N2  =  6.017766  x  102s3 
g2N3  =  7.137658  x  105s2 

g2N4  =  1‘3620729  x  1C)7s 
g2N5  =  4* 7818957  x  108 

82D1  =  1>0s4 

0  3 

82D2  =  9*2  x  10  s 
g2D3  =  2.60192  x  loV 
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g2D4  =  1.5488  x  105s 
g2D5  °  1.024  x  107 


7.  CONCLUSIONS 

The  program  described  in  this  technical  memorandum  for  computing  the 
inverse  Laplace  transforms  and  plotting  their  time  history  response  as  a  function  of 
real  time  has  the  following  merits: 

1.  Uncomplicated  data  loading. 

2.  Produces  equations  in  exact  inverse  Laplace  form. 

3.  Produces  an  accurate  numerical  printout  of  selected  points. 

4.  Produces  a  graphical  display  for  visual  description  of  system  response. 

The  main  objective  of  this  program,  which  is  to  provide  expanded  systems 
analysis  capability,  has  been  realized.  Transient  solutions  of  high-order  differential 
equations  being  driven  by  a  wide  variety  of  forcing  functions  can  be  generated 
rapidly. 
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TABLE  1. 


INPUT  FOR  EXAMPLE  1 
«i(s)  =  \ 


Mirror 

Of  Input 

Columns 

Description 

Input  for  Example 

1. 

(1-80) 

Identification 

0.000 

2.500 

(1-10) (11-20) 

Running  Time 
Start-Stop 

0000 

(1-80) 

Separator  (Blank) 

1 

0 

3.000 

11.000 

(1-5) (6-20) 

Degree-Coefficient 
of  Numerator 

3 

2 

1 

0 

1.000 

8.000 

15 . 000 
0.000 

(1-5) (6-20) 

Degree-Coefficient 
of  Denominator 

5.000 

3.000 

0.000 

0.000 

1 

1 

(1-15) (16-30) 
(31-45) 

Pole  Root 

Real  Part-Imaginary 
Part-Multiplicity 

0000 

(1-80) 

Separator  (Blank) 
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TABLE  2. 


■k 

INPUT  FOR  EXAMPLE  2 


g-^s) 


1 

s 


Mirror  of 

Input 

Columns 

Description 

Input  for  Example  2. 

(1-80) 

Identification 

0.000 

5.000 

(1-10) (11-20) 

Running  Time 

0.000 

1.200 

Start-Stop 

0.000 

.500 

.500 

2.000 

0000 

(1-80) 

Separator  (Blank) 

2 

1. 340E-03 

(1-5) (6-20) 

Degree-Coefficient 

1 

2.067E-01 

of  Numerator 

0 

1.000 

6 

6 . 9981E-07 

(1-5) (6-20) 

Degree-Coefficient 

5 

1. 1552908E- 

-04 

of  Denominator 

4 

1.68E-03 

3 

3. 828E-02 

2 

2.593E-01 

1 

1.000 

0 

0.000 

4.1258865 

4.2939554 

1 

(1-15) (16-30) 

Pole  Root 

2.6595745 

16.89876 

1 

(31-45) 

Real-Par t-Imaginary- 

151.5151 

0.000 

1 

Multiplicity 

0000 

(1-80) 

Separator  (Blank) 

See  example  on  page  165  of  Reference  6. 


Nixon,  Floyd  E. ,  'Principles  of  Automatic  Controls,  Prentice  Hall, 
Inc. ,  New  York,  1954. 
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TABLE  3. 


INPUT  FOR  EXAMPLE  3 

1_ 

2 
s 


gX(s) 


Mirror  Of  Input  Columns 

Input  for  Example  3.  (1-80) 


Description 

Identification 


0.000  .1000  (1-10) (11-20)  Running  Time 

0.000  .0100  Start-Stop 

0.000  10.000 


0000 


(1-80)  Separator  (Blank) 


4  24.6665  (1-5) (6-20) 

3  601.7766 

2  7 . 137658E05 

1  1. 3620729E07 

0  4. 7818957E08 


Degree-Coefficient 
of  Numerator 


6  1.000  (1-5) (6-20)  Degree-Coefficient 

5  9.200  of  Denominator 

4  2.60192E04 

3  1. 5488E05 

2  1.024E07 

1  0.000 

0  0.000 


3.000 

19.77372 

1  (1-15) (16-30) 

Pole  Root 

1.600 

159.9992 

1  (31-45) 

Real  Part-Imaginary 
Part-Multiplicity 

0000 

(1-80) 

Separator  (Blank) 
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APPENDIX  A 


LISTING  OF  SOURCE  DECK 
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TT  19.  .24 


NOV. 14,69  8RL ESC  FORTRAN  2  ANO  4 
♦SA046  J.  ANDRESE  BLDG  390  X2611  TRANSIENT  RESPONSE  2.  4 


C  CANNOT  EXCEEOTHE  FOLLOW  ING. .  .3R0  OROER  POLE  AT  ORIGIN, 2N0  ORDER  POLE 

1 

C  NOT  AT  ORIGIN,  1ST  OROER  COMPLEX  POLE.  ACCEPTS  STABLE  ANO  UNSTABL 

2 

C  SYSTEMS.  ACCEPTS  UP  TO  15  COMPUTING  INTERVALS. 

3 

C 

PUT  8LANK  CARO  AFTER  PRESELECTEO  TIME  AND  AXIS  VALUES  IN  DATA 

4 

C  PUT  BLANK  CARD  AFTER  EACH  PROBLEM  IN  OATA 

5 

C  FINOING  THE  INVERSE  LAPLACE  TRANSFORM  OF  A  QUOTIENT  OF  TWO 

6 

C  POLYNOMIALS  (INVERSE  TRANSFORM  PARAMETERS) 

7 

REAL  PDTOP ( 60  ) 

8 

REAL  C0EF,N0,Nl,01,02,03,C0D(60) , CON (60)  ,DIF 

9 

REAL  RE ( 60 ) , I M ( 60 ) 

10 

REAL  N  ( 5 1 ) , NP ( 50 ) ,0(52) ,0P(51 ) , OPP ( 50 ) , DPPP ( 49 ) 

11 

INTEGER  ON.DNPl, CD.ODPl,  I  ,K 

12 

REAL  PALPHI60) ,P8ETA(60) ,PA0R8(60) ,PTHRC(60) 

13 

INTEGER  KK ( 60  ) 

14 

INTEGER  COMENT( 20  ) 

15 

INTEGER  ITT 

16 

REAL  TST A ( l 5 ) ,TSTOP( 15) ,YSTA( 15) ,YSTOP( 15) 

17 

90 

FORMAT ( 1H  , IP2F18.7, I4,5X,1H( , E18.7,3H)+( , E 1 8. 7 , 3H ) * T  ) 

18 

91 

FORMAT ( 1H  , 1P2E18. 7, I 4, 5 X , IH ( , El  8. 7 , 1 5H ) *EXP ( -ALPH A*T ) ) 

19 

92 

FORMAT  (  1H  ,1P2E18.7,I4,5X,1H(,E18.7,3H)  +  (,E18.7,5H)*TM, 

20 

1E18.7.5H  )*T*T  ) 

21 

93 

FORMAT ( 1H  , 1P2E18.7, I4.5X, IH( , El  8. 7, 17H)*EXP(-ALPHA*T)  +  ( 

,  E  1 8  •  7 , 

22 

1 17H) *T*EXP(-ALPHA*T)  ) 

23 

95 

FORMAT ( 1H  .1P2E18.7, 14, 5X, 

24 

1 1 H ( , E 1 8 . 7 ,30H)*EXP(— ALPHA*T)*SIN( ( 8E T A*T ) ♦ ( , t 1 8. 7 , 2H ) ) ) 

25 

ICSE=1 

26 

999 

CONTINUE 

27 

WRITE (6,  150)  I  CSE 

28 

150 

FORMAT ( 1H1 ,/////////////////////////// ///,60X, 12HPR08LEM 

NUM. , 14 ) 

29 

RE ADI  5 , 888 )  COMENT 

30 

888 

FORMA  T ( 20A4 ) 

31 

WR I TE ( 6 , 77  7  )  COMENT 

32 

777 

FORMAT ( 1H  ,//,20X,20A4) 

33 

99 

JTA=1 

34 

9999 

CONTINUE 

35 

C  READ  REAL  TIME  COMPUTING  INTERVALS 

36 

RE AD( 5  *  101 )  TSTA ( JTA ) , TST0P( JTA) , YSTA (JTA ), Y STOP ( JTA ) 

37 

101 

FORMAT ( 4F 1 0 • 0 ) 

38 

I F ( TSTA ( JTA) .EQ. 0.0. AND. T STOP! JTA). EQ.0.0  ) GOT 09998 

39 

J  T  A* J  T  A+ l 

40 

IF ( JTA* GT. 15)  GO  TO  99 

41 

GO  TO  9999 

42 

9998 

CONTINUE 

43 

C 

CLEAR  PARAMETERS 

44 

DO  500  LCR*lt60 

45 

RE  <  LCR ) =0.0 

46 

CO  N  (  LC  R  )  *=0 . 0 

47 

IM(LCR)=0.0 

48 

COD ( L  CK ) =0 . 0 

49 

PD  TOP ( LCR ) =0 . 0 

50 

PAL  PHI LC  R ) =0 . 0 

51 

KK ( LCR ) =0 • 0 

52 

P AORB ( LCR ) =0 • 0 

53 

P8ET  A ( LCK ) *0 • 0 

54 

500 

P  THRC I LC  R ) *0 . 0 

55 

DO  510  M* 1 » 52 

56 

0 ( M ) *0. 0 

57 

COO ( M ) =0 • 0 

58 

510 

CONTINUE 

59 
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00  511  M= 1 , 5 1 

60 

CON( M ) *0 .0 

61 

RE ( M) =0.0 

62 

IM(M)*0.0 

63 

N(M)=0.0 

64 

DP  I M ) *0 .0 

65 

511 

CONTINUE 

66 

DO  512  M= 1 , 50 

67 

NP ( M ) =0. 0  0 

68 

DPP I M  )  =  0 .0 

69 

312 

CONTINUE 

70 

DO  513  M*  1,49 

71 

OPPP ( M) =0. 0 

72 

513 

CONTINUE 

73 

WRITEI6, 15) 

74 

15 

FORMAT ( 1H1,7X,5H ORDER, 

11X 

♦4HC0EF  f 

/) 

75 

C 

READ  COEF.  OF  HIGHEST 

POWER 

OF  S 

IN  NUMERATOR 

76 

RE AD ( 5 , 1  )  K , COEF 

77 

1 

FORMAT! I 5, F 15 .0) 

78 

WR I TE ( 6 , 1 1 1 ) 

79 

111 

FORMAT ( 1H  , 7X , 18HTERMS 

IN 

NUMERATOR, / ) 

80 

WR  I  TE ( 6 , 19  )  K.COEF 

81 

19 

FORMAT ( 1H  ,5X, 14 , 5  X, 1P2E18. 

7) 

82 

I  11=0 

83 

KNXK  ♦ 1 

84 

I  =5  1-K 

85 

KP  l  =  K  ♦ 1 

86 

CON ( KP  l ) =  COEF 

87 

ON*  I 

88 

ONP 1  =  DN  + l 

89 

N ( I  )  *COEF 

90 

C 

READ  REST  OF  COEF.  OF 

NUMERATOR 

91 

GO  TO  60 

92 

2 

RE AO ( 5 » 1 )  K  » COE F 

93 

WR I TE ( 6 1 19  )  K  f  CO  EF 

94 

1=5 1-K 

95 

KP 1 =K ♦ 1 

96 

CON ( KP 1 ) =COEF 

97 

N (  I  )  =COE F 

98 

C  ZERO  COEF  CHECK, START  COEF 

OF 

DENOMINATOR 

99 

60 

I F ( K.NE • 0 ) GO  TO  2 

100 

C 

READ  COEF  OF  HIGH.  POWER 

OF 

S  IN 

DEN. 

101 

WR I T  E ( 6 1 112) 

102 

112 

FORMAT ( 1 H  ,  /  / , 7X , 2 OH TERMS 

IN  DENOMINATOR,/) 

103 

KE  AO ( 5  t 1 ) K » COEF 

104 

WR  I  TE ( 6  « 19  )  K  t  CO  EF 

105 

KD=K+l 

106 

I*52-K 

107 

KP  l  =  K+ 1 

108 

COL)  (  KP  1  )  =COEF 

109 

C 

DEGREE  OF  DEN. 

110 

DD=  I 

111 

DDPl=DD+l 

112 

D (  I  )  =COE F 

113 

REST  OF  COEF.  OF  DEN. 

114 

3 

READ(5»1)K»C0EF 

1  15 

WR I TE ( 6 « 19)  K  tCOEF 

116 

1 3  52- K 

117 

KP 1*K+  1 

29 

l  18 

COO ( KP 1 ) =  COEF 

119 

o  o 


0(1) *C0EF  120 

IF(K.NE.0)G0  TO  3  121 

END  READ  I NG  OF  COEF, START  DERIVATIVE  CALCULATIONS  122 

DERIVATIVE  OF  OEN.  123 

00  4  I *DD»  5 1  124 

A I *52- I  125 

4  DPI  I )*AI*D(  I  )  126 

C  2ND  DERIVATIVE  OF  DEN.  127 

00  5  I=0D,50  128 

A  I *5 1-  I  129 

5  DPP ( I ) = A  I *UP I  I  )  130 

C  3RD  DER.  OF  DEN.  131 

00  6  I*DD,  49  132 

A  1=50- I  133 

6  DPPPI I ) = A  I *DP  P ( I  )  134 

C  GET  DER.  OF  NUM.  135 

DO  7  I  =  UN»  50  136 

AI=51-I  137 

7  NP ( I ) =A I *N ( I  )  138 

WR I TE ( 6 »  115)  139 

115  FORMAT ( IH  , // //, 0X  ,5HALPH A, 12X ,4 HBET A , IOX , 1HK,  I9X ,  140 

120HINVERSE  LAPLACE  FORM,//)  14  L 

C  FIRST  OROER  POLE  AT  ORIGIN  142 

IF(DI52).EQ. 0.0. AND. 0(51). NE.O.OJGO  TO  LL  143 

GO  TO  10  144 

11  ALPHA=0. 0  1 A5 

BETA=0. 0  146 

K* 1  147 

AORB=NI 51) /D( 51)  140 

TH0RC=0.0  149 

WRITE! 6,90 ) ALPHA .BETA, K, AORB, THORC  150 

I  I  I  *  I  I  I  +  1  151 

PALPH( I  I  I )  =  A  L  P  HA  t  KK(III)=Kt  P AORB ( I  I  I ) = AOR Bt  152 

PBETA ( I  I  I  )  =Bfc  TAJ  PTHRC ( I  1 1 )*THORC  153 

GO  TO  40  154 

10  IF(DI52).EQ.O.O.ANO.DI51).EO.O.O.ANO.D(50).NE.O.O)GOT013  155 

G0TG12  156 

C  SECOND  ORDER  POLE  AT  ORIGIN  157 

13  ALPHA=0. 0  158 

BET  A=0. 0  159 

K=2  160 

AORB= I D ( 50 ) *N I  50 )-D( 49 ) *N ( 5 1 ) ) /D ( 50 ) **2  161 

THORC=N ( 51 ) /D ( 50 )  162 

WR I TE (6, 90) ALPHA , B ET A , K , AOR B, T HOKC  163 

111=111+1  164 

PALPHJ  I  I  I )=ALPHA*  KK(III)*K$  P AORB(  I  I  I ) = AOR B»  165 

PBETAJ II  I)*BETAJ  PTHRC ( I  1 1 ) *THORC  166 

GO  TO  40  167 

12  IF(D(52).EQ.O.O.AND.D(51).EQ.0.O.AND.D(5O).EQ.O.0)GOTO14  160 

C  NO  POLES  AT  ORIGIN  -  GO  TO  READ  POLE  CO.  169 

GO  TO  40  170 

C  THIRD  ORDER  POLE  AT  ORIGIN  171 

14  ALPHA=0 . 0  172 

BETA=0.0  173 

K=3  174 

A0RB=(D(49)*0(49)*N(49)-D(49)*D( 47 ) *N ( 5 1 )-D ( 49 ) *D ( 40 ) *N ( 50 )  175 

1*NC51  )*D(48  ) *D ( 4 8 )  )/ (0(49 )*D( 49) *D( 49) )  176 

IH0RC=(DI49)*N(50)-D(48)*N(51)  ) / ( D ( 49 ) *0 ( 49 )  )  177 

DT0P  =  N(51)/(2.0*D(49)  )  170 

WRITE (6, 92) ALPHA,BETA,K,A0R6,TH0RC,0T0P  179 

30 


III*III+1  lflO 

PALPHI III )*ALPHA$  KK  {  I  1 1  )  =K $PADR8 I 1 1  I ) =ADRB  181 

PBETAI I  I  I )*BETA*  PTHRC ( I  I  I ) *THDKC *  POTOP ( III  )  =  DTOP  1B2 

C  POLES  AT  ORIGIN  TESTED  FOR, READ  OTHER  POLES  1B3 

40  READI5,2D)ALPHA,BETA,K  184 

20  FORMAT ( 2F15.D, 15 )  185 

WRITEI6.9)  ALPHA, BETA, K  186 

9  FORMAT! 1H0, 1P2E1B. 7, I4.2E18.7)  1B7 

I F ( ALPHA. EQ.D.O.  AND. BETA. EO.D.  )  GO  TO  1D9  188 

IF(BETA.NE.D.D)GDTD5D  189 

I F ( K. NE • 2 ) GO  TO  3D  190 

ND=N ( DN )  191 

DD  21  I =DNP 1 ,51  192 

21  ND*-ALPHA*NO+NII )  193 

IFIDNPl.GT .51 )NO=N(DN)  194 

N 1 *NP ( DN )  195 

DO  22  I =DNP 1 , 50  196 

22  N1*-ALPHA*N1+NP(  I  )  197 

IFIDNPl.GT. 5D)N1*NP(DN)  19B 

D2*DPP I DD )  199 

DD  23  I =DDP 1 , 5D  200 

23  D2*-ALPHA*D2  +  DPP  I  I  )  201 

IF!DDP1.GT.5D)D2=DPPIDD)  202 

D3=DPPP I DD )  203 

DD  24  I =DDP 1 , 49  204 

24  D3=-ALPHA*D3+DPPP I  I)  205 

IF (DDP1 .GT .49 )D3=DPPP! DD)  206 

ADRB* I D2/2.D*N1-ND*D3/6.D)/ 102**2/4.0)  20 7 

THDRC=2.D*ND/D2  20B 

WRITE (6,93) ALPHA , BET  A , K , AORB , T HD RC  2  09 

111=111+1  210 

PALPHI I  I  I )=ALPHA*  KK I  I  I  I ) =K$  PAORBI I  I l ) =AORB*  211 

PBETAI I l I )*BETAi  PTHRC I  I  1 1 ) *THORC  212 

GO  TO  40  213 

3D  IFIBETA.NE.D.OIGO  TD  5D  214 

NO=N( DN )  215 

DD  31  I =DNP 1 ,51  216 

31  ND  =  -ALPHA*NO  +  N 1 1  )  217 

IFIDNPl.GT. 51)ND*NIDN)  218 

D 1*DP I DD  )  219 

UO  32  I =DDP 1,51  220 

32  D1=-ALPHA*D1  +  DPI  I)  221 

IF!0DP1.GT.51)D1=DP!DD)  222 

ADRB=ND/D1  223 

THDRC=D. D  224 

WRITE (6,91 ) ALPHA , BETA, K, ADRB, T HO RC  225 

111=111+1  226 

PALPHI I  I  I )*ALPHA  $  KK(III)  =  K  *  P ADRB I  I  I  I )  =  ADR B  227 

PBETAI I  I  I )*BETA  t  PT HRC I  I  1 1 ) *T HO RC  22B 

GD  TD  4D  229 

C  COMPLEX  PDLES  230 

50  TDTNR*CDN! 1 ) i  TDTNI*D.i  TDT DR=CDO I  2 ) i  TDTD I =0. D  231 

A=-ALPHA  t  B  =  BETA  %  AX*l.Dt  AY*O.D  232 

DD  3DD  I R=  1 , KD  233 

REI  IR)=A*AX-B*AY  234 

IMI IR)=B*AX+A*AY  235 

AX*RE I  I R )  236 

A Y*  IM I  I R )  237 

3D0  CONTINUE  238 

DD  301  I N= 1 , KN  239 
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I NP 1= I N*  1  240 

IFIINPl.  GT.KN)G0  TO  305  241 

TOTNR=CON(  IN+1  )*RE  (INH-TOTNR  242 

301  TOTNI-CONl IN* 1 )* IM ( I N ) +TOTN I  243 

305  CONTINUE  244 

DO  302  I D=  l  *  KD  245 

DI F= I D*  l  246 

I DP2  = 10*2  247 

I F (  10P2.GT.KD)G0  TO  306  24B 

T  OTDR=COD ( ID*2)*REIID)*DIF*T0TDR  249 

302  TOTDI-CODI ID*2)*IM( I D > *D I F* TOTDI  250 

306  CONTINUE  251 

AMD=SQRT( TOTDR**2*(TOTDI**2) )  252 

AMN*SQRT (T0TNR**2* IT0TNI**2) )  253 

ANGD=ATAN2( TOTDR.TOTDI )  254 

ANGN=ATAN2( TOTNR.TOTNl  )  255 

A0RB=2.0*AMN/AMD  256 

THORC- A NGN- A NGD  257 

Y1M*A0RB*S1N(TH0RC)  258 

YRE*  AORB*COSI THORC)  259 

TH0KC*ATAN2( YRE.YIM)  260 

C  IF  IN  FIRST  OR  FOURTH  QUAD,  LEAVE  IT  261 

C  OTHER  WISE  ADO  OR  SUB  PI  TO  GET  SUCH,  AND  REVERSE  262 

C  SIGN  TO  COMPENSATE  263 

IFIABSITHORO.LE. 1.57079632)  GO  TO  453  264 

IF(TH0RC)450,451,452  265 

450  TH0RC*TH0RC*3. 14159265  266 

GOTO  451  267 

452  TH0RC=TH0RC-3. 1415926  26B 

451  AORB=-AORB  269 

453  CONTINUE  270 

WRI TE( 6, 95) ALPHA, BETA, K.AORB,  THORC  271 

I  I  1*1  I  1  +  1  272 

PALPHI 1 1  I )  =  ALPHA$  KKI I  I  I  )  =  K$  PAORBI I  I  I )*AORB$  273 

PBETAI I  I  I )=BETA$  PTHRC ( l 1 1 ) * THORC  274 

GO  TO  40  275 

109  CONTINUE  276 

00  1100  I TT= 1 » 15  277 

CALL  PLOT(PALPH,PBETA,KK,PAORB,PTHRC,TSTA,TSTOP,ITT,PDTOP,  278 

1YSTA.YST0P)  279 

IF(TSTA(ITT*1).EQ.0.0.AND.TST0PI ITT*1).EQ.0.0)G0  TO  1200  280 

1100  CONTINUE  281 

1200  CONTINUE  2B2 

ICSE=  ICSE-*-l  283 

GOT  0999  284 

END  285 

C  END  COMPUTATIONS, START  GRAPHICAL  DISPLAY  286 

C  PROGRAM  TO  CHOOSE  TIME  INTERVAL  AND  GRAPH  RESULTS  287 

SUBROUTINE  PLOT t PALPH, PBET A ,KK , P AORB, PTHRC , T ST A, TSTOP ,  2BB 

1ITT,PDT0P,YSTA,YST0P)  2B9 

REAL  TSTAI 15) .TSTOPI 15 ),YSTAI 15) ,YSTOP( 15)  290 

INTEGER  ITT  291 

REAL  VALTI 15)  292 

REAL  FT  t 15  )  , TIME  I  15 )  293 

REAL  PALPH(60),PBETAI60),PA0RB(60),PTHRCI60)  294 

REAL  PDTOP 1 60 )  295 

INTEGER  KKI60) , K ( 60 ) , I ,L,J  296 

REAL  ALPHAI60)  ,BETA(60), A0RBI60) , THORC ( 60  I , E , F , T , TF , FP  297 

REAL  ALPHAS, ALPHAL, BETAS, BETAL .DELTAT , TEND, DELTA  298 

INTEGER  IS  299 
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REAL  GR APH ( 60  » 1 30 )  300 

DATA  BLANK/ 1H  /,  DOT/LH./,  X/LHX/  301 

TEND=TSTOP< ITT)  302 

BP0*-99999999999.  303 

SNE=99999999999.  $ZAP=SNE  304 

A  I *0. S  B J=0  .  S  F I *0. OS  FJ=0.0  305 

DO  430  NCL=1,15  306 

FT (  NC  L  )  =  0 . 0  307 

TIMEtNCL )=0.0  308 

430  VALT!NCL)=0.0  309 

FT(l)*1.0S  FT ( 2 ) -2  *0  S  FT!3)=3.0S  FT!4)=4.0S  FT!5)*5.0  310 

VALT! 1  )  =  TSTA!  ITT )  ill 

GEST=0.0  312 

BEST=0. 0  313 

TEST=0.0  314 

C  ARRAY  FOR  GRAPH  315 

DO  612  1=1,60  316 

DO  612  J= l » 130  317 

612  GRAPH! I , J)=BLANK  31B 

C  DOTS  FOR  AXES  319 

DO  667  1=10,60  320 

667  GK APH ( I » 10 ) =DOT  321 

DO  670  1=10,12  322 

DO  670  J=10, 130, 10  323 

670  GRAPH! I , J)=DOT  324 

DO  66B  J  =  10 , 1 30  325 

66B  GRAPH! 10, J )=DOT  326 

C  WRITE  HEADING  327 

WR I T  E ! 6 »  999B )  32B 

999B  FORMAT!  IH1,  15X.IHT , 13X,3HFOT , 10X, 5HDELT  A, 7X, 6HDELTAT , // )  3  29 

DO  501  L  =  1 , 60  330 

M=L  331 

ALPHA ! M)=PALPH!M )  332 

BETA!M)=PBETA!M)  333 

K! M)=KK! M)  334 

AORB!M)=PAORB!M)  335 

THORC !M)=PTHRC!M )  336 

IFIK(M).EQ.O)  GO  TO  64  337 

501  CONTINUE  338 

64  L=L-1  339 

C  GET  ALPHAS, ALPHAL, IS, BETAS, BET AL  340 

I  S= 1  341 

ALPHAS=ABS! ALPHA !  1  ) )  342 

ALPHAL=ABS ! ALPHA ! 1 ) )  343 

BETAS=ABS!BETA!1  )  )  344 

BETAL  =  ABS! BETA! 1  )  )  345 

DO  23  1=2, L  346 

IFIALPHAS.LT.ABS  IALPHA! I )  )  )G0  TO  24  347 

ALPHAS=ABS ! ALPHA ! I ) )  34B 

IS=I  349 

24  ALPHAL=AMAX1!ALPHAL,ABS!ALPHA! I) ) )  350 

IFIBETA! I ) . EO. 0. 0 ) GO  TO  23  351 

BETAS=AMIN1(BETAS,ABS!BETA! I)) )  352 

BETAL=AMAX1!BETAL,ABS!BETA!I>> )  353 

23  CONTINUE  354 

327  CONTINUE  355 

C  FIND  INITIAL  STEP  SIZE, DEL  TAT  356 

IF! ALPHAL. GT.BETAL  )G0  TO  25  357 

DELTAT=. 1/BETAL  358 

GO  TO  26  359 
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25  DELTAT*. 1/ALPHAL 

26  CONTINUE 

T*  TSTA<  I TT) 

C  EVALUATE  FUNCTION 

67  F*0  •  0 

DO  78  I*1,L 
I F ( BE  T  A (  I  )  )80, 81,80 
C  COMPLEX-NOT  ZERO  BETA 

80  F*F  +  AORB ( I ) *EXP( -ALPHA ( I ) *T ) *S IN  I  BETA  II) ♦T  +  THOKCI  I  ) ) 
GO  TO  78 

C  REAL  -  BETA  ZERO 

81  E« EXP ( -ALPHA (  I  )♦ T) 

F=F+AORB ( I ) *E 

I F ( K ( I ) .EQ.2)F*F+TH0RC( I ) ♦T ♦£ 

I  F  (  K  (  I  )  .EQ.3)F*F+PDT0P(  I  )*T*T*E 
78  CONTINUE 

IF ( YSTA( I TT) .NE.O.O.AND. YSTOPI  I T T ) .NE . 0 . 0  )G0 TO 326 
IF( TEST. EQ. 1.0)  GO  TO  326 
DO  411  NTE*  l ,4 
FT ( NTE ) =  FT ( N  TE  +1) 

TIME(NTE)*TIME(NTE  +  l) 

411  CONTINUE 
FT ( 5 ) *F 

T I  ME ( 5 )  =  T 

IF(T.GT.TSTOP( ITT) )G0  TO  333 
C  DETERMINE  SCALING  OF  AXIS 
SNE-AMIN1 ( SNE,F) 

BP0*AMAX1 (BPO,F) 

DO  412  NTE=2,5 
DJAA=ABS(FT (NTE-l)-FT(NTE ) ) 

IFIDJAA.GT. ( ABS ( FT ( 1 ) )♦. 00000001) )GOTO  49 

412  CONTINUE 

TEND- TIME! l)$  TEST*l.O 
GO  TO  327 

333  TEND  =  TSTOP(  ITTJS  TESTM.O 
GO  TO  327 
326  CONTINUE 

IF(BEST.EQ.l.O)  GO  TO  399 

I F ( YSTA ( I TT ) .EQ.O.O. AND. Y STOP!  ITT  )  .EQ. 0 . 0 ) GO  TO 305 
BPO*  Y  ST  0  P (  ITT) 

SNE- Y  ST  A (  ITT) 

305  CONTINUE 

FI*50.0/(BP0-SNE ) 

FJ*120.0/(  T END-T  S T  A(  ITT)) 

BE  ST*  1 . 0 
399  CONTINUE 
TF*F 

WRITE (6, 72) T, TF, DELTA, DEL  TAT 
72  FORMAT ( l H  , 4X , 2F 1 5 .4 , 6X , 2 F9 .4 ) 

A I  * ( F I ♦ ( TF-SNE ) )  + 1 0.  0 
IF (AI.GT. 100000.  ) GOT  0450 
I  *  I  F  I  X  (  A  I  ) 

IF(GEST.EQ.l.0)G0T039 
ZP  T*BPO 

SZRO* (BP0-SNEJ/50.0 
DO  350  K T* l , 50 
ZAP*AMIN1(ZAP,ABS(ZPT) ) 

ZPT*ZPT-SZRO 
350  CONTINUE 

Z  I  *  (  F  I  *  ( Z  AP-SNE )  )+10.0 
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JA* I F I  X ( Z I )  420 

DO  92  J=10,  130,5  421 

92  GRAPH! JA, J)=OOT  422 

GE  ST* 1  423 

39  CONTINUE  424 

8J*!FJ*IT-TSTA!ITT)) )«-10.0  425 

IFJ8J.GT. 100000. 00)G0T0  450  426 

J=IFIX!BJ)  427 

IF(TF.EO.SNE)  1*10  428 

1 F ( TF.EQ.BPO)  1=60  429 

1F( T.EQ.TSTAI ITT ) )  J*10  430 

IF(T.EO.TENO)  J*130  431 

1 F I  I . GE . 10.AN0.I .LE.60.AN0.J.GE. 10. AND. J .LE. 130 )G0  TO  400  432 

GO  TO  450  433 

400  GRAPH ( I , J ) =X  434 

450  CONTINUE  435 

I  F I T .GE  . TENO ) GO  TO  666  436 

C  FINO  0ERIVAT1VE  AT  T*T  AND  T*T*DELTAT  437 

49  FP*0. 0  438 

00  178  1=1, L  440 

E  =  EXP1-ALPHA1  1  )*T)  441 

IF 18ETA11) ) 180,181,180  442 

C  COMPLEX  -  BETA  NOT  ZERO  443 

180  FP*FP«-A0R8I  1  )*E*COSI  BETA!  I)  *T*THORC!  1  )  )*BETA!I)  444 

1  -A0RB1 1  )*ALPHA!  1  ) *E*SIN(BETA( I) *T*THORC! 1 ) )  445 

GO  TO  17B  446 

C  REAL  -  BETA  ZERO  447 

1B1  FP*FP- ALPHA!  I  ) *A0R8!  I )  *E  448 

I F ( K ( I ) • EQ.2)FP=FP+TH0RC! I ) *E* ! I . O-ALPHA I  I ) *  T )  44  9 

I F ( K ( I  1 .EQ.3)FP=FP*P0T0P! I )*E*!2.0*T-ALPHAI I)*T*T)  450 

178  CONTINUE  451 

AX*2.0*FP  452 

Y=T *DEL  T  AT  453 

FP*0. 0  454 

DO  27B  1*1, L  455 

E*EXP!-ALPHAI 1 )*Y)  456 

I F ! BE  TA !  1 ) 1280,281,280  457 

280  FP=FP*AORB !  I  ) *E*CO  SI  BETA!  I  )*Y«-THORCI  I  )  ) *  BET  A !  I  )  458 

1  -A0R8! I  1  * ALPHA! I)*E*SIN(BETA(I) *Y+TH0RC 1 1) )  459 

GO  TO  27B  460 

2B1  FP*FP-ALPHA! 1  )*AORB!  I )*E  461 

I F ( K ( I ) .EQ.2)FP*FP+TH0RC ! 1 ) *E* ! I .O-ALPHA ! I ) * Y I  462 

1F1KI 1 ).EQ.3)FP=FP+PDT0PI 1 ) *E* I  2 . 0*Y- ALPHA ! I  )*Y*Y)  463 

278  CONTINUE  464 

FP*2 • 0*FP  465 

OELT A=A8 S ! AT AN ! AX ) -AT AN ! FP) )  466 

C  END  DERV.  EVAL.  467 

1FIDELTA.GE.0.051G0  TO  50  468 

T  =  T ♦DEL  T  AT  469 

IF!T.GE.TEND)T*TENO  470 

DELTAT*1. 500*0 EL  TAT  471 

GO  TO  67  472 

50  IF!OELTA.GE.O.  15.AND.DELTAT.GE. ! TSTOP! ITTl-TSTA! ITT) )/  473 

1120. 0)G0  TO  6C  474 

T=T*DELTAT  475 

1F!T.GE.TEN0)T=TEN0  476 

GO  TO  67  477 

60  DELTAT*DELTAT/2.0  478 

GO  TO  49  479 

666  WR1TEI6,120)TSTA!ITT),TEN0,TST0P!1TT)  480 
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120  FORMAT! 1H1 * 11HSTART  T l ME* , F 10. 4, 1  OX, 17HACTUAL  STOP  TIME=,F10.4,  481 

110X, 19HPR0P0SE0  STOP  TIME*,F10.4)  482 

WR 1TE ( 6, 8 )  483 

8  FORMAT ( 1H  ,4X,4HF(T) ,/)  484 

VALY=BPO  485 

SCAY=(BPO-SNE  J/50.0  486 

DO  29  11=1,51  487 

1*61-11  488 

WRITE (6, 30) VALY,  (GRAPH! I , J ) , J* 10 ,  1 30 )  489 

VALY*VALY-SCAY  490 

29  CONTINUE  491 

30  FORMAT ( 1 H  , F8 . 3, l X , 1 21 Al )  492 

SCAX=( TEND-VALT!  1) 1/12.00  493 

DO  680  LAX*2, 13  494 

680  VALT(LAX)=VALT(LAX-1 )*SCAX  495 

WRITE (6, 5) (V ALT (LAX) ,LAX=l, 13)  496 

5  FORMAT ( 1H  , 2X , 12 F 10. 3, F9 . 3 )  497 

WR ITE ( 6 , 7 )  498 

7  FORMAT ( 1H  , / , 50X , 1 4HRE AL  TIME(SEC))  499 

C  TAKE  ANOTHER  SET  OF  DATA  500 

RETURN  501 

END  502 
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APPENDIX  B 


FLOW  CHART  OF  PROGRAM 
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*099^  CONTINUE 

*C  READ  RfcAL  T (^E  COMPUTING  INTERVALS 

♦  REA0(!>t10I)  TSTAI  JTA)  t  TST0PI  JTA),YSTA(  JTA)  ,YSTCP(  JTA) 

«  101  FORMAT (4P10.0) 


IF(T$TA(jrA).Eg.0.0.AN0.rST0P(JTA).EG.0.0  )G0F09998 
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TABLE  5.  EXAMPLE  1  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  2.5  SECONDS)  FOR  FIGURE  2 
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TABLE  6.  EXAMPLE  1  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  3 


< 

i 

a 

x 

(JJ 


o 

I 


CL 


< 

I 

a. 

x 

ai 

* 

o 

o 


oj 

t. 

O 

o 

o 

O 

r 

o 


Q. 

X 

OJ 


a 

i 

UJ 

o 

o 

n 

o 

r> 

o 

o 

• 

U*\ 

I 


X 


o 


o  o  o  o 

o  o  n  n 


o  o  o 

o  r> 


r~  — <  o  o 

o  o  cj  o  n  o 


UJ  UJ 
Ci  o 

ll  a  o 

ut  r  o 

o  a:  n  o 

1 1  /~i  o  o 

j-  o  o 

<r  o 

CC  •  • 

JJ  O  ^ 

i: 

3 


UJ  Lit  UJ  UJ 

o  o  o  r. 

a  oo  o  O 

r  r  c  -  <*  o 

►—  o  o  o  o 

<r  o  a  k  ”i 

7*  Ci  o  o  r> 

—  C  O  O  lA 

r  •  ■  •  • 

O  — l  O  (M  ^ 

CJ 

o 


uj  UJ  uj  uj 

o  o  o  o 

o  o  o  o 

o  o  o  o 

o  o  o  o 

O  o  r  i  C> 

c>o  r>  o 

or  or" 

•  •  mm 

no  on 


UJ  UJ  UJ 

O  o  o 

r>  o  o 

O  c.  o 

no  o 

Cj  o  r. 

o  o  o 

r  n  n 

*  •  • 

on  n 


i  o  n  ">  o  -•  o  n  o  "» 
o  o  o  no  o  o  o  o  o 


re  l/>  on 

my  y 

o  tv  nr 

tv  uj  ^  o  uj  m  n  ^  n 

1  j—  j— 


<r 


CL 


Uj  UJ  UJ 

r,  n  -> 
non 
o  o  o 
o  r>  o 
o  o  n 
o  o  n 
o  o  o 


UJ  LL  UJ 

o  o  n 
o  o  o 
O  O  o 
r»  r>  o 

non 
n  o  o 
o  o  o 


UJ  UJ  UJ  UJ 

r>  n  o  o 

non  o 
o  o  o  o 

o  r  o  o 

o  n  o  n 

o  r>  o  o 

o  o  o  o 


•  ••  •••  •  •  •  • 

moo  ■— <  o  in  ia  o  o 


75 


TABLE  7.  EXAMPLE  1  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  5.0  SECONDS)  FOR  FIGURE  3 


_J  OOUMninuMry(VrjNNNfNNN(NJ(NMMN(\jrgn^Olf»^O^'0O^O<fl’4-'0l/'NJ3 

lu  rg*i4,>44-4’4'n4rgcgrgrgrgcgrgrgrgrgcgcgrgmminn»*-Mr-ininininininxr,-'OX'4’ 

o  r,  O'TooonoooonocooooodonoooH^(\j(\ji\j(\rM(\j^irMrN0' 


fnx<-iO<»ng<l'riX'X>omr'»XX^4,inXr,-OC'-i4,4,rr>X4’XXin<~iO>-'<^^4)4,mr- 
i—  o  n  ->oo.~i~iooo^— <— <—i.-i->ooo.-)'nr*i<n'-')Oor->o“'>r>oooor'>->or'» 


Q 


n  cr  ia 
o  X  — I 
►—  f~'  _r\  m 

"5  O  O  ^ 


j-  A  tO  >t  •*>  <rj 

rg  rg  in  in  — <  4 
N  O'  J-  X  O  — < 
f\I  f\»  rrv  rr>  -4- 


o  x  x  m 
4-  m  O  X 
rvl  ^  4  4 
4-  4-  4-4- 


4-^X4- 
O  <X  4  X 

in  in  in  in 

4-  4-  4-4- 


O  5*  -*  X 
in  m  rg  O 
in  in  in  4- 
4-  4-  4-4- 


n. 

x 

4- 

4- 


»*•  -O  O  O'  (VI  ff> 
h  4  4  4-  <0 

4-  m  rg  o  r*>  rg 
4-  4-  4-  4-  ro  rr> 


o  in  x  r>  o 
hiaoian 
*0  O  O  N  O' 
rg  rg  ^  O 


fA  O'  N  N 

in  —i  x  rg  >o 
n  a  rgn  n 
o  o  o  o  o 


o  O 

r.  rg 

n  o 


.■>  in 
in  at 
o  r, 


r-jooinriinoinoinoin 
O  A^N’ANONlANOfM 
4©miANO'N4-O0CHm 


o  in 
in  r» 
in  n* 

4-  4- 


o  in  ^ 
O  rg  n 
C  rg  4 

in  in  m 


n  m  r>  45 
N  H  A  lA 
\fl  O  Ml  ® 
l A  A  O  O 


40 


in  m  43  rr 
in  X  rg  ® 
r-  4  no  m 
s  n  m  m 


(\J  A  0D  r-l 

in  — •  a-  4- 
p-1  MM  30 
®  o  m  in 


n 

x 

x 

o 


in 

4- 

in 


4- 


-3 

ro 

o 


,-i^^HrgfvjrgpgMi4in 


7#. 


TABLE  8.  EXAMPLE  1  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  4, 
G(s)  =  (1  +  0.5s)/(l  +  0.3s  +  0.25s2) 
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TABLE  9.  EXAMPLE  1  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  5.0  SECONDS)  FOR  FIGURE  4 
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TABLE  9  (Continued) 
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TERMS  IN  NUMERATOR 
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TABLE  11.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  5.0  SECONDS)  FOR  FIGURE  5 
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TABLE  11  (Continued) 
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TABLE  12.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  to  1.2  SECONDS)  FOR  FIGURE  6 
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TABLE  12  (Continued) 
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TABLE  12  (Continued) 
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TABLE  13.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  0.5  SECOND)  FOR  FIGURE  7 
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TABLE  13  (Continued) 
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TABLE  14.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0.5  TO  2.0  SECONDS)  FOR  FIGURE  8 
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TABLE  14  (Continued) 
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TABLE  15.  EXAMPLE  2  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  9 
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TABLE  16.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  3.0  SECONDS)  FOR  FIGURE  9 
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-  oo  o  t  o  o  *»  o  o  o  o  o  o  o  o  o  o  'lnoMHiMiMNNmmmnj’tf  ^xMrxNtrpXixiX!  sC>cjic^h-h-cwxxxxxtro'(T'Os 
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TABLE  16  (Continued) 


i/\f\J(VJNNNfSINJ,-<iAir\NNNMM(NNMO'CD0D®eDCD®fOinNNNNa5rOi^O'O'<T'J-Nn-C 

.-i<MiMiMiM<M<MiMAA<MiMf--i«-t«-««-t,-i«-i,-if-»,-iiMAjrMiMiMiM>3-AAAAA^rr*-r-'r>ooA»trwA 

oooooooaooooooaooaaooooooooooooooooo^^f-Hf-HiMr^irN 


s^cD<tecir.cr^<o^in>j-ir\<o 
nrtHiNjNNmnon  .-i  r>  a 


r-NtO'OOAAarsAAO'AAtMA 
>t  N  fO  ^  (T  OOOCOOr-r^P'-AA 
O'  0'0'X>0AAr5>0iMA»-«r-A<M 
ooor»r>Aor>o>.-i,-i^Ar>A 


O'  AA»-*»-*iA'©«-«r>^'GO<*'iMO' 
A  O'  O'  f 

or>oooo»-i0f-ior>'*)0  "> 


A  nJ- 

r^*  r- 

A  A 


GO 

r- 

A 


IM 

oo 

p- 

A 


S  «0  N  h  ^  O 

^  S  N  ^  O'  O  H 

N  «0  in  >t  N  (T  O' 

fO  'O  fO  (SJ 


O'  A  »-* 

N  N  CO 

r-  r-  a 
im  im  ru 


^^fns-0O'0njf0OSf0ONuM)N^^in>fON'0 

<t^0'S'fl-0UMn>fNNH-nm0'ON-0N®^f00'O 

'OAAAAAAAAA^-'tAtM^r'O'GOaO^-n-'O^-A 

NNNNNMfSJNNfSINWNNnjNMHHHHHH^ 


A 

r- 

o 


^  p-  r> 
O'  O 
CO  h-  lA 
o  r>  r> 


O' 

IM 

A 

O' 


.-h  o  <r  r- 

00  ^  (A  «A 
NON  ^ 

O'  o  n  r. 


lA  ,3-  A  —t 
O'  n  in  jo 
-A  O'  •— I  a 

n~)  (A  *—4  t—4 


O'  P-  A  nJ-  A  «— < 
A  A  O'  IM  A  00 
IM  A  -3-  >0  N  CC 
IM  f\l  f\l  <M  IM  IM 


O'  r-  a  O'  ao 

A  A  O'  00  f- 

iH  IM  m  lA  GO 

fO  (A  A  A  A 


A  A 

r—*  >3 

>3 


AAAO'^J'AA.-iA'J- 

>*-AtMA®^AAX>f- 

■J-AlAl/NAAA^-r'-f'- 


NO1® 

1  A  A  A  A  CO  A 
A  A  s3-  A  A  IM 
0O  O'  O  ^  IM  ^  A 


O 

O 


•— <  >H  fH  iH  >H  H  iH  tH  H  H  H  H  H  M  rH  H  H  H  H  M  M  H  H  H  ^  (M  (M  IM  IM  IM  <A 
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TABLE  17.  EXAMPLE  2  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  10 


X 

* 

a 


cc 

LU 

> 


< 

I 

Q. 


* 


r\i 

o 


o 

I 


m 

in 

m 

r- 

<0 


2 

* 


* 

< 

x 

a. 


CO 

I 

LU 

00 

o 

3D 

m 

r- 

co 

00 

• 

fM 

I 


o 

o 


r- 

in 

*0 

r- 

o 

O' 


CD 

O 

I 

LL> 

o 

o 

m 


o 


in 


of 


m  — •  r» 
o  r>  o 

i  i 

LL)  UJ  LU 

o  <~i  r, 
-~>  o  o 
r-  r>  r. 
e.  o  o 
o  r-  o 
O  c» 
moo 


p-  4  m  fM 
o  o  o  o 
till 
LU  tU  UJ  LU 

or sj  oo  o 

O  «0  P*M 

o  od  n  in 
^  m  h-  in 
on  O'  in 
O'  n  -*  r*» 
O'  ^  n  u 


o  — <  — •  4 


-g  o  o  o 
o  o  o  o 

I 

UJ  UJ  UJ  LU 

f-i  n  r»  o 
o  o  o  o 
o  on  o  o 
4  oo  m  o 
in  n  O'  n 
r-  in  in  cj 
m  m  m  f' 


“Yl  H  f\J  U 


O' 

UJ 

o 

r> 


ru 


o 


<o  m  j-  f<i  ru  u  o 


f\i  fM 

r-<  ^ 

— 

o 

o  o 

o  o 

f-<  ^ 

o  o 

o 

o  o 

O  O 

c.  o 

Cj  o 

o 

UJ  UJ 

UJ  UJ 

UJ  UJ 

UJ  UJ 

LU 

o  o 

sO  >0 

1  ^ 

o  o 

o 

O  O 

^  o 

O'  O' 

o  o 

o 

r>  o 

fM  fM 

o 

o  n 

o 

< 

o  o 

in  in 

O  O 

o 

L— 

o  o 

m  fn 

O'  O' 

c  o 

c- 

UJ 

o  o 

O'  O' 

o  o 

o  o 

o 

cc 

r  o 

fM  fM 

>0  4> 

o  r 

o 

•  a 

•  • 

•  • 

•  • 

• 

r,  r> 

4  4 

^-1  —• 

O  “> 

o 

* ) 

O  *  > 

fM  fM  O 

o 

n  o 

o  o 

o  o 

O  O  C' 

o 

UJ  UJ 

UJ  u.» 

UJ  LL 

LU  LL  LL 

UJ 

o  o 

O'  O' 

o  o 

o 

o  o 

fM  fM 

m  m 

m  m  o 

o 

O  O 

co  oo 

*0  nO 

^ 

w 

o  r> 

p*»  p- 

O  O 

fM  fM  O 

c 

o  o 

in  <n 

JD  *0 

m  m  o 

o 

o  o 

fM  fM 

in  in 

—  o 

o 

o  o 

r-t  <—• 

'O  'O 

in  in  o 

o 

•  • 

•  • 

•  • 

•  •  • 

• 

rM  H 

4  >t 

fM  fM 

H  H  O 

o 
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TABLE  18.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  3.0  SECONDS)  FOR  FIGURE  10 


*}■ 

't 


LLI 

O 


o  m  rg  m 

'3^^rgriinP-*-<-Oininpgrgpgpgpg 

OOOOOOOOOOOOOOOO 


r^r>-f-.p-oinxpgrg^F^^,^r\irvgr 

pgpgpgpgO'XPJ.t^rgrgpgpgxpg® 

OOOOOOOOCOOOOOO^ 


<4’P'>p-p-r^P~r~0*^^P'P-m®pg'*-*- 
inpgpgrgpgrgrg.T'CC7vrr'moo>o<7'<7' 
OOOOOOOC-CO-^^Ps/rn^  «£«£> 


pg  r  oc  ^  m  M  ^  c  Cr 
Lnc>ooco-<'j->r>r 
<  '0r3OC>OOO-3^'T 

I-  O  <*5  o  J  O  O  O  O  o 


UJ 

a 


•-ip-fninxcino(7'^«'Oin 
m^srvjCT'HcOin^OO'^'O 
KiCT'O^noaD^'O^ri^ri 
»-iO’-*’-*~4^0OOC00  0 


C^MVf  C'C-O^'fO'^S 

^^rsjQT'ircO^H^xJ’^O' 

oxxotf'P-Pr'o-<®r^CP:7' 

o  ->  ^  o  *5  n  o  o  o  3  o  o 


p'lOP-r^x^rp-x.-gp-rg'-  rg 

sOPgor-a'-^0'OP-«^rom,^ 

^  -o 

000^3^0^ -*>*'  O  '"> 


,  v»  O  ) 

o  o  o  o  o 
o  o  o  o  o 
o  o  o  o  o 


O  •-)  z>  * 

o  o  o  o  o 

r>  o  o  o  o 

o  o  o  o  o 


O  <0  CP  Cj  O 

pg  m  in  ^  m 

o  o  o  o  -< 

o  o  o  o  o 


^  pg  in  x  ^ 

to  <  O'  O' 

^  pg  m  m  < 

O  O  O  O  O 


•*■  in  p»  m 
cr  -o  •*  cp  -* 

OO^^N 


(3  in  ^  ^  iTi 

«o  <o  m  r»  y> 

m  in  r-  od  a* 

pg  pg  pg  »N  pg 


pg  33  oo  pg  rg 

p-  m  ^  cp  in 

^  ^  ^  P 

(*)  (V)  PO  CP 


soo>o«fin^g'^^s  ^  c.  cr  m 
(ppgg’inming-no^o'oonw 
p-xxxx®xxr^,orocrrPin 


o 

r- 

>0 

p— 4 

>r 

p- 

pg 

in 

X 

n 

pg 

O' 

n 

pg 

O' 

«o 

m 

o 

p- 

pg 

o 

m 

n 

<0 

p- 

X 

O 

pg 

'T 

m 

>0 

p- 

X 

O' 

o 

pH 

p- 

p- 

pH 

m 

* 

pg 

X 

p- 

r~j 

o 

„> 

m 

n 

X 

■n 

pg 

O' 

* 

r- 

o 

pg 

in 

X 

o 

m 

<o 

(7* 

H 

o 

<? 

pg 

•o 

X 

pg 

-o 

m 

p> 

X 

m 

pg 

O' 

■o 

* 

p«H 

pg 

*■ 

pH 

X 

* 

pg 

m 

O 

o 

o 

C 

o 

o 

pg 

m 

* 

p- 

X 

o 

H 

pg 

fO 

m 

-o 

p- 

X 

o 

pg 

* 

<7* 

m 

X 

pg 

X 

X 

<0 

pp 

<T 

p- 

O' 

pg 

in 

X 

pg 

X 

p- 

'T 

in 

-o 

pg 

r 

►- 

o 

c 

o 

o 

o 

o 

o 

"V 

o 

o 

o 

^■4 

H 

w-i 

^-4 

H 

pg 

pg 

pg 

pg 

m 

tr i 

*■ 

* 

X 

in 

<o 

<o 

p- 

r- 

p- 

r- 

X 

X 

X 

O' 

CT' 

o 

pg 

pn 

m 

X 

m 

o 

—i  ^  pg  m 


94 


TABLE  19.  EXAMPLE  2  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  11 
G(s)  =  1/(0. Ills2  +  0.266s  +  1) 


** 

rH  H 

o 

o 

r> 

o  o 

o  c 

o 

Cl 

o 

o  o 

o  o 

o 

o 

o 

UJ 

UJ 

Ul  UJ 

Ul  UJ 

UJ 

UJ 

UJ 

H  O 

03 

>7 

o 

rg 

r-4 

O 

CD 

® 

>o  «o 

H  «-H 

o 

o 

o 

0  0  0 

o 

o 

O 

o 

o  o 

r>  o 

O 

O' 

sO  'O 

O'  O' 

O 

D 

o 

1  1 

1 

i 

1 

i 

1 

i 

i 

gi 

«o 

rg  rg 

c  o 

o 

r> 

.  ) 

LL>  U  LLf 

UJ 

UJ 

UJ 

UJ 

Ul 

UJ 

Uj  UJ 

nD 

vC 

in  m 

l-H  l-H 

o 

o 

6 

o  r»  o 

H 

m 

r- 

X 

in 

s7 

o  o 

F— 

O' 

O' 

m  m 

O'  O' 

o 

o 

L 

LL 

o  o  o 

ex: 

O' 

CD 

h- 

m 

r- 

<7 

o  o 

LU 

>7 

s7 

O'  O' 

o  o 

o 

o 

o 

UJ 

o  r  o 

ri 

an 

CD 

■J- 

rv, 

o 

■J- 

r>  o 

X 

r- 

r- 

rg  rg 

gD  o 

n 

c 

<r 

o 

ce 

o  o  o 

F— 

03 

l-H 

in 

u- 

m 

rg 

sO  o 

• 

• 

•  • 

•  • 

« 

• 

• 

1  ) 

.*  • 

or*-  i 

<r 

<7 

OJ 

iH 

n 

O 

m 

o  m 

rg 

rg 

>7  n7 

<— <  <— < 

rj 

L~ 

F— 

^  oo 

7 

r- 

c 

X 

>o 

X 

m  o 

< 

O  o  o 

H 

r- 

rr 

iH 

® 

O 

iH 

rg  O 

ce 

It  •  • 

a: 

• 

• 

• 

• 

• 

• 

o  • 

Ul 

<"*  rg  <"* 

r- 

rg 

s7 

rg 

n 

ac 

7 

"3 

UJ 

n 

O 

O  J) 

O  o 

rg 

rg 

73 

7 

O 

o 

O 

o  o 

o  o 

O  O  O 

O 

7" 

UJ 

UJ 

UJ  UJ 

UJ  UJ 

UJ 

UJ 

Ul 

UJ 

M 

M 

o 

j 

O  O' 

«o 

v7 

o 

O 

< 

o 

o 

rg  rg 

m  m 

m  m 

o 

o 

re 

l/> 

X 

O' 

O' 

CD  ® 

sO  o 

o 

o 

UJ 

X 

y 

CL 

o 

o 

r-  r- 

g3  no 

rg  rg  o 

o 

o 

re 

re 

-J 

o 

<r 

m  m 

'O  nO 

m  m 

o 

o 

re 

UJ 

NH  O 

UJ 

r- 

X 

>7 

o 

rg 

o 

< 

o 

o 

rg  rg 

m  m 

i-h 

i-h 

o 

o 

o 

F— 

F— 

rg 

rg 

i-h  i-h 

nO  nO 

m 

m 

o 

• 

• 

•  • 

•  • 

• 

• 

• 

• 

>— « 

>7  '7 

rg  rg 

o 

o 
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TABLE  20.  EXAMPLE  2  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  3.0  SECONDS)  FOR  FIGURE  11 


UJ 

o 


r*-  ii  rsj 
O  *  TM 
O  O  O  O 
o  o  o  o 


OOOOOOOOOOOOOOOOO^OCjOOOOOOOGOOOOOOO 


•0 

o 


rsj  m 

*  -h 

lA  00 

o  c 


^<\i>rr-  oo  <*  r-  f-  f- 

rsjacf»-m>opnm— 

^Hf\jHOOO<'Cr 


^OChNSiTO^NMh 

<1  (T  J  O  >  U  O  ^  W  (M  -m  N  IT 

F—  ^.i  1  o  ^  O  O  C  h  N  N  N 


>omm>oo^mma'f\i'0^-aop*-4A 
-*  ~>  ~i  o  o  «*•.  o  o  o  t  *-«  ~i  (\j  >r 


ON'CNNU'C5'NirNH(?OX 
0-Oriel-MrHOH(?>0(fiO'oO 
^<03'ONl/>lTiriA'T<OiAW^ff' 
IT  N1  H  H  ,  )  O  •■■>  o  O  ^  O  O  O  O  *"5 


r^p-LACT'^-xiP-^^HO'r'rrg^r- 
'10>rH>rsO'(7'3rg>0'M^onr 
Hf'JW3'(\J7'(NJP-^nHrvi  >T  lA  O 


ui 

o 


ooooooooC'(Mir^r.iffitMO>o^'J'3Q^^'J'(N'Nju:M^x^r\jr-(a'ir(T.r^^Mir-a'Hj->rir>T^'TNr-ii>ifNjoo-'>T^N 

oc'5ooooooco^fN/miAr^O'Ooofs)0'iA'r(N^iA«ooA4A^fri^-<<joor-Lr^oooc'r\)crfn®<T'mof^>x<^(7'o'^o'*1 

ooo^ov.»'‘ooo,'iooc,,:'O^-^<NiAa''0<7'C'”wO'^oooooo^'O^a'ac0Df»-^^Lrc)Fr»'TiAAiALf'ALr^b'C^0'0 


^^rgrjfNjrgfNjrgfNjfsjfNjrgrg^^^^^^^^^ 


!  I  I  (  I  I  I  I  I  I  I  I  I 


TKOh^S  _ _  „  ,  .  _ 

5  ^  ^  m  CNlHfjffU/'JN 

O  O  O  O  O  'J^fNjrr 

O  “  o  O  C  o  ■;  o 


fM^oo^(^>o(7'(Nj<'p''OC'^^x  THfsjfyjnrn<j,^(r<o-c^s®»(?ij^N'fN>o(r>0(v)irxT'-'N4,i/'N  ®  o 

“ - -  -  ~  S\  rf'yj  ^T\  1(MX>rO'£lfNjQ04'O^lNm^  ‘j  »  >a  c  O  *■*  •*« 


(DHNlvOO'lM'TN.O  — <  ^  X)  jH(M(\Mirri^^ir>o^)p.NCOIIO'0^rvJ^r-^(r>ONir 

W  IT  «0  N  C  '.•>  H  (NJ  cr  ^  ff  yj  (T,  'UMX^C>0(M®4'O^N®^  j^'>^7(f''f<0t?0tri7' 

c*  c»  f  ;o^-M»MHH(M(Njffi^<<<r«Tir^(niri/Mr'0'0<o«C'0'O^P'Xocr. 


r,  cnr-<7'-~'r\j>rAr'-C?‘^ 

^  IT  -/l  /l  IT  'A  O 
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TABLE  20  (Continued) 


h-r-r-h-r-r-r-r- 

OOOCIOOOO 


r-  r- 

\f\  in 
pg  pg 
O  O 


iO'flCMflPl-O®!  CO  CO  >4  >4 
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— <  — <  ^  ^  o  d  o  o  o  a  n  -*  o  ^  o  d  o  «*>  o  jo'ionooD^ 
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.fl-OCNSNN® 


^  i/>  -4  r-  o 

t  pg  o  in  -o  o 
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TABLE  21.  EXAMPLE  3  -  INPUT  AND  ANALYTIC  TIME  SOLUTION  FOR  FIGURE  12 

G(s)  -  ]/s2 
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TABLE  22.  EXAMPLE  3  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  0.1  SECOND)  FOR  FIGURE  12 
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TABLE  23.  EXAMPLE  3  -  TABULATION  OF  RESPONSE  VERSUS  TIME 
(0  TO  0.01  SECOND)  FOR  FIGURE  13 
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START  TIKE*  *0000  ACTUAL  STCP  T(Mc*  2.5000  PRCP0SE0  S TCP  TIME*  2.5000 
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FIGURE  2.  TRANSIENT  RESPONSE  OF  EXAMPLE  1  TO  A  STEP  INPUT 
t  «  0  TO  2.5  SECONDS 
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FIGURE  3.  TRANSIENT  RESPONSE  OF  EXAMPLE  1  TO  AN  EXPONENTIAL  LAG, 
t  =  0  TO  5.0  SECONDS 
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FIGURE  4.  TRANSIENT  RESPONSE  OF  EXAMPLE  1  TO  A  DAMPED  SINUSOID 
t  »  0  TO  5.0  SECONDS 
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FIGURE  5.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  STEP  INPUT 
t  -  0  TO  5.0  SECONDS 
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FIGURE  6.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  STEP  INPUT 
t  =  0  TO  1.2  SECONDS 
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FIGURE  7.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  STEP  INPUT 
t  -  0  TO  0.5  SECOND 
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FIGURE  8.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  STEP  INPUT 
t  =  0.5  TO  2.0  SECONDS 
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FIGURE  9.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  AN  EXPONENTIAL  LAG, 
t  *  0  TO  3.0  SECONDS 
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FIGURE  10.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  SECOND  ORDER 
EXPONENTIAL  LAG,  t  =  0  TO  3.0  SECONDS 
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FIGURE  11.  TRANSIENT  RESPONSE  OF  EXAMPLE  2  TO  A  DAMPED  SINUSOID 
t  =  0  TO  3.0  SECONDS 
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FIGURE  12.  TRANSIENT  RESPONSE  OF  EXAMPLE  3  TO  A  RAMP  INPUT, 

t  *  0  TO  0.1  SECOND 
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FIGURE  13.  TRANSIENT  RESPONSE  OF  EXAMPLE  3  TO  A  RAMP  INPUT. 

t  =  0  TO  0.01  SECOND 
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FIGURE  14.  TRANSIENT  RESPONSE  OF  EXAMPLE  3  TO  A  RAMP  INPUT. 
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